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ABSTRACT 



We perform a detailed investigation of sources from the Cold Cores Catalogue of Planck Objects (C3P0). Our goal is to probe the reliability 
of the detections, validate the separation between warm and cold dust emission components, provide the first glimpse at the nature, internal 
morphology and physical characterictics of the P/awc^-detected sources. We focus on a sub-sample of ten sources from the C3PO list, selected to 
sample different environments, from high latitude cirrus to nearby (150pc) and remote (2kpc) molecular complexes. We present Planck surface 
brightness maps and derive the dust temperature, emissivity spectral index, and column densities of the fields. With the help of higher resolution 
Herschel and AKARI continuum observations and molecular line data, we investigate the morphology of the sources and the properties of the 
substructures at scales below the Planck beam size. The cold clumps detected by Planck are found to be located on large-scale filamentary (or 
H cometary) structures that extend up to 20 pc in the remote sources. The thickness of these filaments ranges between 0.3 and 3pc, for column 
H y densities A^H2 ~ 0.1 to 1.6 x 10^^ cm~^, and with linear mass density covering a broad range, between 15 and 400 M© pc~^ The dust temperatures 
r* ■ are low (between 10 and 15 K) and the Planck cold clumps correspond to local minima of the line-of-sight averaged dust temperature in these 
^ ^ ■ fields. These low temperatures are confirmed when AKARI and Herschel data are added to the spectral energy distributions. Herschel data reveal 
a wealth of substructure within the Planck cold clumps. In all cases (except two sources harbouring young stellar objects), the substructures are 
found to be colder, with temperatures as low as 7 K. Molecular line observations provide gas column densities which are consistent with those 
^ inferred from the dust. The linewidths are all supra- thermal, providing large virial linear mass densities in the range 10 to 300 M© pc~^ , comparable 
- - within factors of a few, to the gas linear mass densities. The analysis of this small set of cold clumps already probes a broad variety of structures 
in the C3P0 sample, probably associated with different evolutionary stages, from cold and starless clumps, to young protostellar objects still 
embedded in their cold surrounding cloud. Because of the all- sky coverage and its sensitivity, Planck is able to detect and locate the coldest spots 
in massive elongated structures that may be the long-searched for progenitors of stellar clusters. 

Key words. ISM: clouds - Infrared: ISM - Submillimeter: ISM - dust, extinction - Stars: formation - Stars: protostars 



1. Introduction 

The main difficulty in understanding star formation lies in the 
vast range of scales involved in the process, including not only 
the gravitationally unstable densest and coldest structures, but 
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also their large-scale environment. The characteristics of star 
formation, described by the stellar mass distribution, the forma- 
tion efficiency, the evolutionary timescales, and by the modes 
of formation (clustered or isolated, spontaneous or triggered) 
are linked to the properties of the cold cores and those of their 
Galactic surroundings. Theory predicts that the initial mass func- 
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tion (IMF) is largely determined during the pre- stellar core 
fragmentation phase, but the latter depends on the properties 
of the pre- stellar cores, particularly the nature of their support 
against self-gravity, and their density and temperature distribu- 
tions. Understanding star formation therefore also requires an 
understanding of the formation and evolution of dense cores. 
Turbulence, magnetic fields, and gravity all contribute at sub- 
parsec scales as well as Galactic scales (see reviews of Bergin 
& Tafalla 2007; McKee & Ostriker 2007; Ward-Thompson et al. 
2007; Crutcher et al. 2009). 

So far, ground-based (e.g., Motte et al. 1998; Johnstone et al. 
2000) and Herschel (e.g., Andre et al. 2010; di Francesco et al. 
2010; Konyves et al. 2010) observations of nearby star-forming 
regions have revealed a core mass spectrum with a slope sim- 
ilar to that of the IMF. This intriguing result calls for broader 
investigations, in particular towards the inner Galaxy, where the 
Galactic environment diff'ers from that of the Solar neighbour- 
hood, and the star formation efficiency is anticipated to be dif- 
ferent. 

The physical properties of cold cores are still poorly known, 
possibly because of the short duration of this phase, but mainly 
because of observational difficulties. The thermal emission of 
cold dust, intrinsically weak and blended with that of warmer 
components, must be sought for in the submillimetre range. A 
combination of continuum and molecular line studies is also 
needed to develop a global view of the clouds, from large scales 
(tens of parsecs) down to the scale of dense cores. During the 
past decade, the development of sensitive continuum and het- 
erodyne detectors at millimetre and submillimetre wavelengths 
has significantly increased our knowledge of the properties of 
the cold dark clouds (see Section 1 in Planck Collaboration 
2011s, hereafter Paper I). The Planck^ and Herschel satellites 
now provide a unique opportunity for studying Galactic dust 
emission and in particular its dense and cold component. The 
first Herschel surveys have already provided new insights into 
cold cores, although limited to specific fields (Andre et al. 2010; 
Bontemps et al. 2010; Konyves et al. 2010; Men'shchikov et al. 
2010; Molinari et al. 2010; Motte et al. 2010; Peretto et al. 2010; 
Schneider et al. 2010; Stutz et al. 2010; Ward-Thompson et al. 
2010;Zavagno et al. 2010). 

The Planck satellite (Tauber et al. 2010; Planck 
Collaboration 2011a) provides complementarity to Herschel 
by carrying out an all- sky survey that is well suited for the 
systematic detection of cold cores. In Paper I, we present the 
first statistical results of the cores from this survey. Combining 
Planck and IRAS 100 jim data, we have built a preliminary cata- 
logue of 10,783 cold cores (the Cold Core Catalogue of Planck 
Objects, C3P0). A sub-sample of the most reliable detections is 
delivered as part of the Planck Early Release Compact Source 
Catalogue (ERCSC, see Planck Collaboration 2011c), i.e. the 
Early Cold Cores (ECC) catalogue. The cores from C3P0 cover 
a wide range in properties, with: temperature, from 7 K to 16 K, 
with a peak around 13 K; density from 30 cm"^ to 10^ cm"^, 
with an average value of 2 x 10^ cm"^; mass from 0.3 to 
2.5 X lO'^M©; and size ranging from 0.2 to 18 pc. The sources are 
found to be extended, and more importantly are elongated, with 
a distribution of aspect ratios extending up to 4 and peaking at 



^ Planck (http://www.esa.int/Planck) is a project of the European 
Space Agency (ESA) with instruments provided by two scientific con- 
sortia funded by ESA member states (in particular the lead countries 
France and Italy), with contributions from NASA (USA) and telescope 
reflectors provided by a collaboration between ESA and a scientific con- 
sortium led and funded by Denmark. 



twice the beam size of the Planck-HFl instrument. As discussed 
in Paper I, these properties better match the definition of 
"clumps" rather than "cores" (see Williams et al. 2000; Bergin 
& Tafalla 2007). Clumps may contain dense cores, which are 
likely the precursors of individual or multiple stars. Thus, the 
main population seen with Planck does not correspond to single 
pre-stellar cores, but ensembles of cold substructures. A detailed 
study of the Planck cold objects requires observations at higher 
resolution. This is the main objective of the Herschel Open 
Time Key Program "Galactic Cold Cores" that is dedicated to 
following up of a sample of Planck cores. 

In this paper we present the first detailed analysis of ten rep- 
resentative sources from the C3P0 catalogue. The targets have 
been selected to cover different types of clump, spanning a wide 
range of mass, temperature, and density, and located in different 
environments, ranging from high latitude cirrus to nearby and 
distant molecular complexes. We combine the Planck data with 
available ancillary data (in particular IRIS, 2MASS din&AKARI) 
and dedicated follow-up observations with Herschel and ground- 
based radio telescopes. Using these ancillary data, we seek to 
confirm the Planck detections and to demonstrate the reliabil- 
ity of the method which used only Planck and IRAS data to 
determine the source parameters that enter the C3P0 catalogue 
(Paper I) and ECC catalogue (Planck Collaboration 2011c,v). 
The higher resolution continuum and line data make it possible 
to examine the internal structure of the Planck sources. They 
provide the first hints about the process of internal fragmenta- 
tion and the physical state of the compact cores at scales below 
the size of the Planck beam. These properties are essential for 
the interpretation of the full C3P0 catalogue and will be a ma- 
jor topic for the projects carrying out follow-up studies of C3P0 
and ECC catalogue sources. 

After describing the observational data set and the source 
sample selection (Section 2), we explain the analysis methods 
used to derive physical properties for the cores (Section 3). The 
main results are presented in Section 4. By fitting spectral enery 
distributions (SEDs) we derive temperatures, emissivity spectral 
indices, and column densities for the cores and the surround- 
ing fields, and estimate the linear mass densities, masses and 
bolometric luminosities of the cores (Section 4.1). The small- 
scale structures within the Planck clumps are studied with the 
help of Herschel and AKARI data (Section 4.3.1), and the gas 
properties of a few cores are derived from molecular line data 
(Section 4.3.3). Following a discussion of the physical charac- 
teristics of the Planck cores, we present our summary and per- 
spectives for the future in Section 6. 

2. Observations 

2.1. The sources 

We have selected from the Cold Core Catalogue of Planck 
Objects (C3P0) ten sources with high reliability, signal-to-noise 
ratio S/N > 8, and low colour temperature, T < 14 K (see 
Paper I). An initial Monte Carlo sampling of the full C3P0 cata- 
logue was performed to prepare a candidate list that covered the 
full range in Galactic position, temperature, ffux, and column 
density. Further selection was made by examining the Planck 
data and ancillary information from IRIS (Miville-Deschenes & 
Lagache 2005), 2MASS extinction (Skrutskie et al. 2006), and 
CO line data from Dame et al. (2001) and NANTEN surveys 
(e.g., Fukui et al. 1999), as weU as IRAS mid AKARI point source 
catalogues and by cross-checking the sources with the SIMBAD 
database. The selected sources represent different large-scale 
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morphologies and environments, including filaments, isolated 
and clustered structures and high-latitude cirrus clouds. One key 
criterion was the knowledge of the source distance, derived ei- 
ther with an extinction method (Marshall et al. 2006), by as- 
sociation with a known molecular cloud complex, or through a 
kinematic distance estimate (see paper I). Nine out of the ten 
sources have already been observed with Herschel as part of 
the Open Time Key Program "Galactic cold cores." In partic- 
ular, the sample includes the three targets observed during the 
Herschel Science Demonstration Phase (SDP, see Juvela et al. 
2010, 201 1). Half of the sources in the sample are in the Planck 
ECC. This is mainly due to the criteria used to select the most 
reliable detections from the C3P0 full catalogue (see Paper I): 
with T < 14 K for the colour temperature corresponding to the 
aperture photometry SED, and a high signal to noise ratio for the 
source detection (S /N > 15). 

The targets are listed in Table 1 and presented in detail in 
the Appendix. The ten sources that we focus on will hereafter be 
referred to by the labels SI through SIO. 

2.2. Planck data 

Planck (Tauber et al. 2010; Planck Collaboration 2011a) is the 
third generation space mission to measure the anisotropy of the 
cosmic microwave background (CMB). It observes the sky in 
nine frequency bands covering 30-857 GHz with high sensitiv- 
ity and angular resolution from 3V to 5'. The Low Frequency 
Instrument (LEI; Mandolesi et al. 2010; Bersanelli et al. 2010; 
Mennella et al. 2011) covers the 30, 44, and 70 GHz bands 
with amplifiers cooled to 20 K. The High Frequency Instrument 
(HEI; Lamarre et al. 2010; Planck HEI Core Team 201 la) covers 
the 100, 143, 217, 353, 545, and 857 GHz bands with bolome- 
ters cooled to O.IK. Polarisation is measured in all but the 
highest two bands (Leahy et al. 2010; Rosset et al. 2010). A 
combination of radiative cooling and three mechanical cool- 
ers produces the temperatures needed for the detectors and op- 
tics (Planck Collaboration 2011b). Two data processing centres 
(DPCs) check and calibrate the data and make maps of the sky 
(Planck HEI Core Team 2011b; Zacchei et al. 2011). Planck's 
sensitivity, angular resolution, and frequency coverage make it a 
powerful instrument for Galactic and extragalactic astrophysics 
as well as cosmology. Early astrophysics results are given in 
Planck Collaboration, 201 Ih-z. 

We use data from the Planck-HFl bands at 857 GHz, 
545 GHz and 353 GHz that cover the main peak of the cold 
dust emission. By restricting ourselves to these highest frequen- 
cies, we can perform the analysis at the best angular resolu- 
tion provided by Planck, i.e., ~ 4.5' full width at half-maximum 
(EWHM). 

Eig. 1 displays the 857 GHz surface brightness maps for each 
1° X 1° field. The cores are located at the centres of the maps and 
show dust emission that ranges from 10 to 220MJy sr"^ They 
are embedded in extended structures whose shape varies from 
large filaments (S4, S5, S6, and S7) to more isolated and ap- 
parently compact morphologies. Most cores are extended and 
elongated compared to the Planck beam. As discussed in Paper 
I, the ellipticity and extension of the cores are not biased by the 
local beam shape. As an illustration, Eig. A. 5 shows the compar- 
ison between the local Point Spread Function provided by the 
FEBeCoP tool (Mitra et al. 2010) with the elliptical Gaussian fit 
of the sub- sample detections. Typical sizes are from 0.2 to 1 1 pc. 
Some cores (S2, S8, and S9) are located near bright, warmer re- 
gions, where the intensities are higher by up to a factor of 10. In 



a few cases, the environment exhibits sharp edges, at large scale 
(see in particular SI, S3 and S7). 

The full set of HEI maps (1° x 1° fields from 857 to 100 GHz) 
for each source are shown in Figs. A.l A. 2 and A. 3. These 
maps have been derived from the HEALPix sky maps (Gorski 
et al. 2005). With the exception of source S2, the cold source 
at the centre of the maps is usually visible down to 143 GHz. 
At 100 GHz the sources become difficult to detect because of the 
falling intensity of the dust spectrum. Only the Musca filament is 
clearly visible even at 100 GHz, although the surface brightness 
of the source is low (~ 0.15 MJy sr"^). 

Comparison with the lOOyum maps from IRIS (Miville- 
Deschenes & Lagache 2005) confirms the conclusion of pre- 
vious submillimetre surveys (using PRONAOS, ARCHEOPS, 
BLAST, and ground-based telescopes) that the cold dust emis- 
sion is not traced by the lOOyum data, but must be studied us- 
ing longer wavelengths. The IRIS maps are dominated by the 
warmer and more extended structure around the cores. This is the 
basis of the source detection method CoCoCoDeT (Cold Core 
Colour Detection Tool), that was described in detail in Montier 
et al. (2010) and in Paper I (Section 2.4). It uses as a template the 
spectrum of the warm emission component estimated from the 
IRIS lOOyum map. Subtraction of the warm component results 
in residual maps of the cold emission component from which 
source fluxes are derived. The detection process is applied inde- 
pendently to the 857, 545, and 353 GHz maps, after smoothing 
HEI and IRIS data to the same angular resolution of 4.5'. The 
main steps of the method and the core extraction process are 
illustrated in the Appendix (in Fig. A.4) for SI. The residual 
lOOyum signal that remains when a background model is sub- 
tracted from the IRIS map (see Paper I) is still needed to con- 
strain the temperature of the cold core SEDs, which typically 
peak in the range 200-300 yum. The analysis of Section 4 will 
also test the accuracy of this procedure. The coordinates, dis- 
tances, fluxes densities, and sizes of the selected Planck cores 
are given in Table 1 . The values listed there are taken directly 
from the C3P0 catalogue (see Paper I and Section 3.1). 

2.3. Herschel observations 

The Herschel photometric observations were carried out with 
the PACS and SPIRE instruments (Pilbratt et al. 2010; Poglitsch 
et al. 2010; Griffin et al. 2010). Three fields, S8, S9, SIO (corre- 
sponding to the source names PCC249, PCC288, and PCC550) 
were observed in November and December 2009 as part of the 
Herschel SDP. The other fields were observed between July and 
September 2010. Most observations were performed separately 
with PACS (100 and 160 yum) and with SPIRE (250, 350 and 
500 yum) in scan mapping mode. Because of the larger field size, 
S5 and S9 were observed in parallel mode using both instru- 
ments simultaneously. The observations employed two orthog- 
onal scanning directions, except for the PACS observations of 
PCC288 where three scanning directions were used. The data 
were reduced with Herschel Interactive Processing Environment 
(HIPE), using the official pipeline with the addition of spe- 
cialised reduction routines to take advantage of the orthogonal 
scans for deglitching PACS data and to remove SPIRE scan base- 
lines. The PACS maps were created using the MADmap algo- 
rithm (Cantalupo et al. 2010). The SPIRE maps are the product 
of direct projection onto the sky and averaging of the time or- 
dered data, with a baseline correction. 

As for most bolometer observations without an absolute cal- 
ibrator, the zero level (or offset) of the PACS and SPIRE data are 
arbitrary. We therefore compared the Herschel and PACS data 
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with the predictions of a model constrained by the Planck and 
IRAS data. The model uses the all-sky dust temperature maps 
described in Planck Collaboration (201 lo) to infer the average 
radiation field intensity for each pixel at the common resolu- 
tion level for the Planck and IRAS data. The DUSTEM model 
(Compiegne et al. 2010) with the above value for the radiation 
field intensity was used to predict the expected brightness in 
the HerschelSVW^ and PACS bands, using the nearest avail- 
able Planck or IRAS band for normalisation and taking into ac- 
count the appropriate colour correction in the Herschel filters. 
The predicted brightness was correlated with the observed maps 
smoothed to the Planck and IRAS resolution over the region 
observed with Herschel and the gain and off'sets were derived 
from this correlation. We have used these gain and off'set inter- 
calibration values in order to convert the maps from Jybeam"^ 
and Jypix"^ (SPIRE and PACS, respectively) into brightness 
units (MJy sr"^). 

2.4. AKARI observations 

The AKARI satelhte (Murakami et al. 2007) has conducted all- 
sky surveys at infrared wavelengths centred at 9yum, 18yum, 
65 //m, 90yum, 140 yum, and 160yum. We use the observations 
made by the FIS instrument in the wide far-infrared bands of 
90yum and 140yum. The accuracy of the calibration is currently 
estimated to be 26% at 90 jim and 33% at 140 yum, and the beam 
sizes of these two bands are ~39'' and 58'', respectively. For 
details of the AKARI far-infrared all sky survey, see Doi et al. 
(2009). 

2.5. Molecular line data 

We have carried out "fast" observations of diff'erent CO iso- 
topic lines in some 60 Planck cold core candidate fields for the 
Herschel follow-up programme, five of which are included in the 
present sample. 

The Onsala 20-m telescope was used for ^^CO /=1^0, 
^^CO /=1^0, and C^^CO /=1^0 observations in December 
2009 and April 2010. An area of a few arcmin in diameter was 
mapped around the position of the Planck source S3 and the 
sources in the SDP fields SB and S9. The Onsala beam size is 
approximately -33'' and the typical rms noise in the ^^CO and 
C^^O spectra was below 0. 1 K per a channel of 0.07 km s"^ . 

The APEX observations of field SIO were made in July 2010. 
The size of the ^^CO /=2^1 map is ~5' and the typical rms 
noise is below 0.2 K. The APEX beam size at 220 GHz is ~ 28". 

The IRAM-30m observations were performed in July and 
October 2010. The EMIR receivers E090 (3 mm) and E230 



(1 mm) were used in parallel with the high resolution correla- 
tor (VESPA). The most important parameters for the diff'erent 
settings are given in the Appendix in Table ??. Maps of 3'x3' 
were performed using on-the-fly (OTF) mode combined with 
frequency switching. A summary of the observations and the 
main results obtained are presented in Section 4.3.3. 

3. Methods 

3.1. Photometry and SEDs 

3.1 .1 . Planck and IRAS photometry 

The method for estimating the photometry for the C3P0 cata- 
logue has been described in detail in Montier et al. (2010) and 
Paper I (Section 2.4). We recall here only the main steps of the 
detection and flux extraction process: 

1 . For each pixel, and for each frequency, the warm background 
colour (Cbkg) is estimated as the median value of the ratio of 
the Planck to 100 yum emission maps (/y//ioo) within a 15' 
radius disc; 

2. For each Planck frequency, the contribution of the warm 
component is obtained by extrapolation from lOOyum 
through = Cbkg x/ioo; 

3. The cold residual map is computed by subtracting the warm 
component from the Planck map; 

4. The cold source detection is performed using a thresholding 
method applied on the cold residual map, with the criterion 
S/N > 4; 

5. The source shape is estimated by fitting a 2D elliptical 
Gaussian to the colour map /gsv/Zioo; 

6. The flux density at lOOyum is derived by fitting an elliptical 
Gaussian plus a polynomial surface for the background; 

7. The warm template at lOOyum is corrected by removing the 
source contribution that was estimated in the previous step; 

8. Aperture photometry at 857, 545, and 353 GHz is performed 
on the cold residual maps, with the aperture determined by 
the source shape from step 5. 

The photometric uncertainties associated with this method 
have been estimated with a Monte Carlo analysis (see Paper I, 
Section 2.5); they are 40% for IRAS 100 yum and 8% in the 
Planck bands. The additional calibration uncertainties to be 
taken into account are 13.5% and 7%, respectively, for IRIS 
(Miville-Deschenes & Lagache 2005) and the HFI bands 857, 
545, and 353 GHz (Planck HFI Core Team 201 lb). As described 
in Paper I, the elliptical Gaussian fit performed in step 5 leads to 
an estimate of the minor and major axis lengths (crum and (XMaj, 
respectively), related to the FWHM values, i.e., ^Min and ^Maj of 
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Fig. 1. P/a/tc^-HFI brightness emission maps at 857 GHz. The colour scale is in MJy sr ^ . Note the large dynamic range of the 
Planck cold clumps from 10 MJy sr'^ (S2) to 250 MJy sr'^ (S9). 



a Gaussian by cr = ^/ V8 ln(2). The fluxes and FWHM values 
are given for each source in Table 1 . 

The 353 GHz band includes some contribution from the CO 
/=3^2 molecular line. The magnitude of this efl'ect has been es- 
timated in Planck HFI Core Team (201 lb) using data from ^^CO 
/=1^0 surveys, together with an estimated average line ratio of 



/=3^2 to /=1^0, and knowledge of the 353 GHz band spec- 
tral response. The derived correction factor is 171 yuK in thermo- 
dynamic temperature for a CO /=1^0 line area of 1 Kkms"^. 
The correction is small, but the exact efl'ect is hard to estimate 
for our sources because of the lack of high spatial resolution 
CO data and because the CO excitation towards cold cores may 
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differ significantly from the average values assumed in deriving 
the above factor. We therefore present results without making a 
correction for the CO contribution, but comment on its possible 
influence later. 



quadratic sum of the noise map and the calibration uncertainties. 
The averaged column densities of the clumps are then derived 
from the observed flux at vq = 857 GHz and the dust colour 
temperature inferred from the fit: 



3.1 .2. Photometry with ancillary data 

When using Herschel and AKARI data, the fluxes are estimated 
at a difl'erent resolution and with direct aperture photometry. 
The maps are convolved to a resolution of 37'' (for Herschel 
and AKARI 90yum), 58'' (for Herschel and AKARI 90yum and 
140yum), or 4.5' (when including Planck and IRIS data). In 
each case the aperture has a radius equal to the FWHM of 
the smoothed data and the background is subtracted using the 
30% quantile value within a reference annulus that extends from 
1.3xFWHMto 1.8XFWHM. In ordernotto be afl'ected by possi- 
ble emission from transiently heated grains, the fits employ only 
data at wavelengths longer than lOOyum. The statistical uncer- 
tainty of the flux values is derived from the surface brightness 
fluctuations in the reference annulus. As above, the calibration 
uncertainty is taken to be 7% for the Planck channels and 13.5% 
for IRIS 100 Jim. For Herschel we assume a 15% calibration un- 
certainty for PACS and 12% for SPIRE. 

3.2. Estimation of temperatures, spectral indices, and optical 
depths 

Since the dust thermal emission is optically thin in the submil- 
limetre range, the source SEDs can be modelled as a modified 
blackbody of the form 



(3) 



TyBy(Tc)O^C\ 



(1) 



where 5 y is the flux density integrated over the clump solid an- 
gle Qci (common to all frequencies), Ty = A^H2y^^H^vo(^/^o)^ is 
the dust optical depth, and By is the Planck function at the dust 
colour temperature, T^. Here is the H2 gas column density, 
= 2.33 the mean mass per particle, mu the mass of the proton, 
KyQ the mass absorption coefficient at frequency vq and JS the dust 
emissivity spectral index. The SED may be rewritten as 

Sy=ABy(T^)y^ (2) 

to separate the three quantities (assumed independent) which can 
be derived from the fit: amplitude A oc Nw^Ky^^; spectral index P\ 
and effective temperature T^. The;^^ minimisation search oper- 
ates in the (A, Tc, P) space and takes into account the colour 
corrections for the Planck bands described in Planck HFI Core 
Team (201 la). In the fitting procedure, we consider only the ffux 
error bars associated with the photometry method, i.e., 40% for 
IRIS and 8% for HFI bands. In a second step, we compute the 
final error bars on Tc and (3 by adding the contribution from 
calibration uncertainties; these have to be considered separately 
because the HFI calibration errors are not independent of ea- 
chother. In Paper I the size of this effect was estimated using a 
grid of Tc and p values, and was found to be small (< 2.5% for Tc 
and < 2% for P) compared to the uncertainties of the flux extrac- 
tion method. The final uncertainties in and p are obtained as a 
quadratic sum of the two contributions. They do not include the 
correlations between and p inherent in the fitting procedure 
itself. 

Applying the same fitting method pixel by pixel to the IRIS 
100 yum and HFI 857, 545, and 353 GHz surface brightness im- 
ages (1° X 1°), we calculate maps of the dust colour tempera- 
ture and the emissivity index. Here the flux error bars are the 



with Ilcl = TrCTMajCTMin. 

The value of Ky is a main source of uncertainty. Large vari- 
ations exist between dust models, depending on properties of 
the dust (see Beckwith et al. 1990; Henning et al. 1995): com- 
position (with or without ice mantles); structure (compact or 
fluffy aggregates); and size. Both models and observations show 
that Ky increases from the diffuse medium to dense and cold re- 
gions by a factor of 3-4 (Ossenkopf & Henning 1994; Kruegel 
& Siebenmorgen 1994; Stepnik et al. 2003; Juvela et al. 2011; 
Planck Collaboration 201 lu). For this study, we adopt the dust 
absorption coefficient of Beckwith et al. (1990) inferred for high- 
density environments (Preibisch et al. 1993; Henning et al. 1995; 
Motte et al. 1998): 



Ky = O.lcm^g-^ (y/lOOOGHz)^, 



(4) 



with a fixed emissivity index of p = 2. The reference frequency 
vo = 857 GHz chosen to estimate is such that the extrapo- 
lation from 1000 GHz with p = 2 (instead of the fitted spectral 
index) introduces an uncertainty on A^H2 of less than 30%. At a 
distance d, the mass of the clump is simply M = Qci J^mud^ 
or 



M: 



Sy,d' 



(5) 



4. Results 



4.1. Physical parameters derived from Planck and IRIS data 

The cold clump SEDs determined using the Planck HFI and IRIS 
flux densities are shown in Fig. 2. As discussed in Section 3.1.1, 
these SEDs are those of the cold dust residual. The values ob- 
tained for the dust temperature and emissivity spectral index 
are shown in the frames of the figure. Using the method de- 
scribed in Section 3.2, we infer the gas column density, bolo- 
metric luminosity, and mass of each clump (Table 2). All these 
parameters (T, p, A^Hz) are values which are averaged over the 
line-of-sight and the extent of the cold clumps, which reach the 
parsec- scale. Denser and/or colder structure is to be expected on 
smaller scales. These averaged temperatures vary between 10.3 
and 14.7 K, and the dust emissivity spectral indices range be- 
tween 1.8 and 2.5. 

The column densities are distributed around a mean of 
7.5 X 10^^ cm"^, with a few of the largest values being above 
~10^^cm"^. These relatively moderate column densities also 
correspond to averaged values over the extent of the clumps. 
The clumps are likely to be heterogeneous, with denser substruc- 
tures that will be studied with higher angular resolution data in 
Section 4.3. 1 . The masses of the sample cover a large range from 
3.5 to 1800 Mq. The highest-mass object, S6, is also the coldest 
source of the sample, and it is quite extended (~2.9 x 6.7 pc); 
these properties make this source an interesting candidate for be- 
ing a high mass star-forming precursor in its early stages. More 
generally, the sizes of the source sample are rather large, in the 
range 0.2-1 1 parsecs, which confirms that they should be better 
classified as "clumps" rather than "cores," according to the ter- 
minology used for nearby molecular clouds (e.g., Williams et al. 
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2000; Motte et al. 2007). Scales typically considered for "dense 
cores" and "clumps" are ~ 0.1 pc and ~ 1 pc, respectively. 

As noted in Paper I, most of the cold clumps found with the 
CoCoCoDet procedure are elongated, with aspect ratios signif- 
icantly larger than unity. The maps of Fig. 1 show that these 
clumps are not isolated, but are, in many cases, substructures of 
long filaments, not always straight. Without direct estimates of 
the density it is difficult to assess that they are actual filaments, 
denser than their environment. However, their brightness con- 
trast above the background is a few, in the case of S4, S5, S6 and 
SIO, suggesting actual density enhancements given their small 
thickness. Their length, estimated from their angular size in the 
maps of Fig. 1 and their distance, ranges from about 5 pc for SIO 
to ISpcfor S5. 

Because of this, we choose to focus not on the clump mass, 
that depends on the size of the clump derived from the detec- 
tion procedure, but on the mass per unit length of the filament at 
the position of the cold clump. This mass per unit length (or 
linear mass density) is given by m = 16.7 Mopc"^ A^2i^Min, 
where A^2i is the clump column density expressed in units of 
10^^ cm"^ and ^Min, the half-power thickness of the clump iden- 
tified in the Planck maps, is taken as the width of the filament. 
The linear mass densities are given in Table 2 for all clumps, in- 
cluding those with an aspect ratio only slightly larger than unity. 
However, they should be considered only as an approximate 
guide to the linear mass density at the location of the clump. 

The linear mass densities vary by a factor ~ 30 and compare 
well with those inferred from molecular lines. The range of val- 
ues (15 to 4OOM0pc"^) is characteristic of regions of massive 
or intermediate-mass star formation, such as Orion and p Oph 
(Hily-Blant et al. 2004). The linear mass densities of the most 
tenuous filaments (< 1 Mopc"^ Falgarone et al. 2001) are not 
found in the present sub-sample of Planck cold clumps. 

Many of the cold clumps detected by Planck therefore ap- 
pear as cold substructures within larger scale filamentary struc- 
tures that have length up to ~ 20 pc in the Planck maps and with 
parsec-scale thickness (S5, S6, S7). The thinnest filaments found 
are only a few parsec long and with thickness of a few 0.1 pc 
(SIO). 

At this stage of the analysis, it is difficult to assess the exact 
nature of the sources, i.e., discriminating between protostellar 
objects and starless clumps. We can however use the tempera- 
ture and bolometric luminosity estimates as first indicators of 
the evolutionary stage of the sources. Such an approach was pro- 
posed by Netterfield et al. (2009) and Roy et al. (2010) in their 
analysis of the cold cores detected with the balloon-borne ex- 
periment, BLAST, in the Cygnus-X and Vela surveys. Among 
the values obtained for the ratio L/M in our sample (ranging 
from 0.13 to 0.91 Lq/M©), it is interesting to note that the three 
highest (with L/M > O.6L0/MO) are associated with T > 14K 
dust sources S2, S8 and S9). As we will see in Section 4.3.1, the 
Planck clumps in these fields harbour warm dust sources, and 
two of them are associated with young stellar objects. 

4.2. Colour temperature and spectral index maps 

Fig. 3 shows the 1° x 1° maps for dust colour temperature 
and dust spectral index that were estimated using the IRAS 
lOOyum data and the three highest frequency HFI channels (see 
Section 3.2) without background removal. The three Herschel 
SDP fields (S8, S9 and SIO) are not presented here, because the 
analysis of T and p was already described in detail in Juvela et al. 
(2011). 
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Fig. 3. Maps of dust colour temperature, emissivity spectral in- 
dex and column density for each clump field (the three SDP 
fields are shown in Juvela et al. (2011) using Herschel maps). 
The A^H2 contours correspond to the ticks on the colour bar. For 
T and the contours are shown respectively at 2K intervals and 
for 13= 1.85. 
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Fig. 2. SEDs and fit parameters (T,/3) obtained combining Planck HEI and IRIS flux densities integrated over the clump in the cold 
residual maps. 



Except for the faint (and not so cold) source S2, a clear 
signature of cold dust emission is visible directly on these tem- 
perature maps at the location of the clumps. The colour temper- 
ature characterizing the line-of- sight integrated intensity in the 
direction of the clump is typically about 14 K, warmer by a few 
Kelvin than the clump temperature deduced from its SED after 
warm background subtraction. The clumps are embedded within 
extended structures in the temperature range 16-18 K. For ref- 
erence, dust in the difl'use medium of the local Galactic neigh- 
bourhood has a temperature of around 17.5 K (Boulanger et al. 
1996; Schlegel et al. 1998). In several cases (S4, S5, S6 and S7) 
the cold clumps appear as the coldest spot (or one of the cold- 
est spots) along filamentary structures, which are already colder 
than the larger scale environment. The temperature and column 
density maps are correlated, with the coldest structures corre- 
sponding, as expected, to the most opaque parts of the field at 
the angular resolution probed by Planck. The sharp edges seen 
in the intensity maps of the fields SI, S3 and S7 are also clearly 
associated with column density sharp transition. This may be the 
signature of shock compression and triggering core formation 
that should be investigated in further studies. 



We have compared the results on column densities with ex- 
tinction measurements. Extinction is calculated with the NICER 
method (Lombardi & Alves 2001) using stars from the 2MASS 
catalogue (Skrutskie et al. 2006). The Ay maps of the fields 
are derived assuming an extinction curve that corresponds to 
Rv=3.l. Preliminary extinction maps are created using nearby 
fields where IRAS 100 jim data indicates a low dust column den- 
sity. The Ay maps obtained with this method are shown in the 
Appendix (in Fig. ??) at 4.5' resolution. A clear signature is vis- 
ible, with an increase by about a factor of 2 in Ay toward the 
clumps. The Ay values derived (~ 2-5 mag) are consistent with 
the previous A^H2 estimates of several times lO^^cm"^, derived 
from the SEDs. 

In most fields the colour temperature and the spectral index 
are seen to be anticorrelated, with high spectral indices being 
found at the location of the temperature minima. The parame- 
ters T and J3 are known to be connected so that any noise present 
in the observations tends to produce a similar anticorrelation. 
However, the spatial coherence of the T and JS maps strongly 
suggests that the results of Fig. 3 are not caused by statistical 
noise. There is still the possibility of systematic errors, but the 
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Source 
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^Min 


^Maj 




m 


M 
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L/M 




(K) 




(pc) 


(pc) 


(cm-2) 


(Mopc-i) 


(Mo) 


(L0) 


(Lo/Mo) 




11 6 + 08 
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4.4 X 10^^ 




D .J 


1 


98 




1 a Q _i_ 1 1 
U.O It 1 . i 


Z.J It u. / 


1 9 


1 s 


O.U X iU 


ID 


1 J.U 


1 1 n 

1 i .u 


n 79 

u. / z 


S3 


11.1+0.7 


2.4 + 0.6 


0.4 


0.5 


1.1 X 10^2 


77 


22.0 


4.3 


0.19 


S4 


11.5 + 0.8 


2.3 + 0.7 


0.7 


1.0 


6.3 X 10^1 


74 


44.0 


11.0 


0.24 


S5 


13.3 + 1.1 


2.0 + 0.6 


2.7 


11.0 


1.7 X 10^1 


77 


520.0 


270.0 


0.52 


S6 


10.3 + 0.6 


2.5 + 0.7 


2.9 


6.7 


8.4 X 10^1 


400 


1800.0 


230.0 


0.13 


S7 


12.5 + 1.0 


2.1+0.6 


0.6 


1.5 


4.9 X 10^1 


49 


45.0 


17.0 


0.37 


S8 


14.0+1.5 


1.8 + 0.7 


1.2 


2.2 


8.1 X 10^1 


160 


210.0 


140.0 


0.65 


S9 


14.7 + 1.4 


1.9 + 0.8 


1.0 


1.9 


1.6 X 10^2 


270 


330.0 


300.0 


0.91 


SIO 


11.2 + 0.7 


2.5 + 0.7 


0.3 


1.0 


1.3x1022 


64 


40.0 


8.4 


0.21 



^ Minor and major axis FWHM values of the Gaussian ellipse, converted to length units. 
^ Linear mass density. 



similarity of the results obtained in separate fields makes this un- 
likely. A pure calibration error would change the J3 and T values 
in a systematic way but would not significantly aff'ect the disper- 
sion of the spectral index values (Juvela et al. 201 1). The line-of- 
sight temperature variations tend to decrease the observed spec- 
tral indices. The eff'ect should be stronger toward the cold clump. 
The fact that ft at those positions is larger than in the diff'use ar- 
eas suggests that the emissivity index of the grains has increased 
in the cores more than what is visible in the maps. Therefore, 
the results support the idea of spatial variation of grain proper- 
ties that could result from temperature-dependent processes in 
the dust emission mechanism (Meny et al. 2007; Boudet et al. 
2005; Coupeaud et al. 201 1, see Paper I discussion). 

The CO correction to 353 GHz band fluxes would typically 
increase the jS values by ~ 10% or less, and the efl'ect on the 
derived colour temperatures is at most a few 0.1 K. The main 
efl'ect is a change in the general amplitude of the values without 
major changes in the morphology of the parameter maps. 

4.3. Analysis using ancillary data 
4.3.1 . Herschel and AKARI maps 

Eor all sources we have higher resolution dust continuum data 
in the form of Herschel (resolution 37'' or better) mid/or AKARI 
(58'' or better) maps. These are used: to examine internal struc- 
ture of the sources; to derive physical parameters of the Planck 
sources and compare them with the results of Section 4. 1 ; and to 
study the properties of the structures found at scales below the 
Planck beam size. 

Eig. 4 shows Herschel 250yum data (resolution ~ 18") for 
all, but one, sources. S2 has not been observed with Herschel 
and we show the AKARI 90yum map instead (~ 39"). Eor com- 
parison, the Planck 857 GHz brightness contours are overplot- 
ted and ellipses show the cold clump sizes as derived from el- 
liptical Gaussian fits on the cold residual maps (cf. Paper I and 
Section 3. 1 . 1). A blow-up of the Herschel 250 yum maps is shown 
in Eig. 5. 

In most fields, a significant amount of structure at scales be- 
low the Planck resolution is visible. The Musca field (SIO) is an 
exception. The Herschel observations there resolve the filament 
but do not reveal further structure within the P/a^c/:-detected 
clump. Interestingly, this field is among the closest to the Sun. 
The other filamentary structures (S4, S5, S6, and S7) break up 
into numerous smaller bright knots, the brightest often coincid- 
ing with the position of the Planck detection. The two cometary 



shaped clouds (SI and S3) similarly harbour a number of bright 
smaller knots and narrow filaments that show up with a large 
contrast (> 2) above an extended background. 



4.3.2. Photometry with Planck, Herschel and AKARI 

We have performed fits on the SEDs of selected substructures 
within the Planck cold clumps to further characterise their inner 
structure. The aperture photometry is performed after subtrac- 
tion of a background level estimated as described in Section 3. 
The radius of the largest aperture is 4.5' (i.e., a diameter of 
540"). Eor the substructures, we use apertures of either 74" 
or 116" diameter. The locations of these apertures are drawn 
in Eig. 5 and listed in Table 3. The SED fits involve data at 
wavelengths longer than lOOyum. This removes the problem of 
a possible contribution from transiently heated grains to the 
lOOyum flux densities, but also reduces the efl'ects of the line-of- 
sight temperature variations. Eor the large aperture SEDs, where 
Planck fluxes are used, the 353 GHz data are corrected for the 
CO emission, when available. This correction has little efl'ect on 
the results of the aperture photometry. 

Eig. 6 illustrates the various SED estimates for the SI and 
S5 fields. It displays the Planck SED as given in the C3P0 cat- 
alogue (therefore corrected for the warm background contribu- 
tion) and the aperture photometry using smaller apertures on the 
brightest substructure of the SI field (Eig. 5). Eor the large aper- 
tures, the flux densities from Planck and Herschel are in good 
agreement with each other; in fact it appears that the assumed 
sources of uncertainty may overestimate the errors. The temper- 
atures obtained for the fixed apertures (Table 3) are the same 
within the error bars as those derived from the Planck and IRIS 
data only (Table 2). It is therefore most encouraging that the es- 
timates based only on Planck and IRIS data do not significantly 
change when the five Herschel bands are added that better cover 
the short wavelength side of the SED maxima. The modified 
blackbody parameters derived from the fits are listed in Table 3. 

One of the most interesting results is the low temperatures 
obtained for the SEDs in the 74" apertures. Except for sources 
S8 and S9, they are all significantly smaller than those obtained 
in the larger apertures, the lowest values being at about 7 K. 
Therefore, in spite of the large telescope beam size, the sensivity 
of Planck is such that the cold clump detection method is actu- 
ally able to pick-up and locate (within a few arcmin) structures 
of very cold dust that are significantly smaller than the beam. 
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Fig. 4. Planck 857 GHz surface brightness contours on Herschel SPIRE maps at 250 yum. The source S2 has not been observed with 
Herschel and the displayed image corresponds to the AKARI 90 jim wide filter. The dashed ellipses correspond to the estimated size 
of the Planck cold clump. 
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Fig. 5. Blow-up of Herschel SPIRE 250 jim maps at the location of the Planck cold clumps localised by their elliptic boundary. 
The different circles and the annuli refer to the apertures adopted for the photometry used in the SEDs. The large aperture has a 9' 
diameter, and the smaller ones correspond to either 74'' (S3 and S5), 116'' (SI, S4, S6, S7, S8 and S9), or 360" (SIO). 
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Table 3. Results from SED fits based on flux densities from aperture photometry (for A > lOOyum data). 



Object 


/ 


b 


Aperture 


Aperture size 






A^H2 




M 




(deg) 


(deg) 


(arcsec) 


(pc) 


(K) 




(cm-2) 


(cm ^) 


(M0) 


Sl-1 


126.61 


24.55 


540 


0.39 


12.59 + 0.61 


2.32 + 0.14 


9.3 X 10 






Sl-1 


126.61 


24.55 


116 


0.08 


10.01 + 1.32 


2.62 + 0.54 


1.9 X 10^^ 






Sl-1 


126.61 


24.55 


74 


0.05 


9.75 + 1.98 


2 67 + 92 


3.0 X 10^^ 


2.7 X 10^ 


0.13 


S2-1 


316.56 


20.68 


540 


1.44 


15.08 + 0.93 


1.67 + 0.17 


1.9 X 10 






S3-2 


131.74 


9.70 


540 


0.52 


11.26 + 0.35 


2.36 + 0.11 


3.3x10^1 






S3-3 


131.74 


9.68 


116 


0.11 


7.09 + 0.74 


3.96 + 0.78 


c a 1 A 22 
5.3 X 10 






S3-4 
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In sources SB and S9, the substructures are found to be 
warmer, with: 19-21 K. As presented in Juvela et al. (2010), 
the PACS lOOyum maps have revealed the presence of bright 
compact sources within the P/a^c/:-detected clumps. These ob- 
jects are likely to be very young stellar objects, still embedded in 
a more extended dense and cold cloud, whose emission is domi- 
nating the submillimetre wavelength range studied with Planck. 
As described in the Appendix, SB and S9 are both located close 
to active star-forming regions, with SB having nearby HII re- 
gions, and OB stellar associations. 

We note that neither SB or S9 has any clear counterpart 
in the 100 yum IRAS map (see Eig. A.3); the emission of the 
young stellar objects is not visible, even although this is the 
same wavelength as for PACS, probably because the sources are 
very faint, and diluted with the other components in the IRAS 
beam. However, for SB, both the 170 [im ISO Serendipity Survey 
(Stickel et al. 2007) and the AKARI EIS survey at 140 yum (Doi 
et al. 2009) showed a faint feature at this position. Shorter wave- 
length data with higher angular resolution would be interesting 
to combine for these particular sources. Among our sample, the 
field SB has been observed by Spitzer, with the 24yum MIPS 
data revealing a number of compact sources. Within the Planck 
clump, the brightest submillimetre peak (PCC2BB-A in Juvela 
et al. 2010) is seen to contain at least four distinct mid-infrared 
sources. These data will be analysed in detail in a future publi- 
cation. 

We have also estimated the column densities within each 
aperture, and the average densities and masses obtained in the 
smallest aperture. In most cases, the derived column densities 
are higher than those estimated from the Planck clump SEDs. 
Moreover, as we found with the linear mass densities estimated 
with Planck data, they vary over a factor of about 30. These val- 
ues should probably be considered as lower limits, because, as 
shown in Eig. 5, the annulus associated with the small apertures 



are close to the source and often include a fraction of signal as- 
sociated with the source itself. We have estimated the impact of 
the background level on the derived column density by using the 
average brightness in the largest annulus instead of the smallest. 
The derived column densities are increased slightly, by up to a 
factor of 2. 

4.3.3. Analysis of molecular line data 

All the Planck clumps observed in ^^CO and C^^O lines have a 
clear line detection, and the small maps show peaks that coincide 
with the cold substructures (Eigs. 7 and B). 

The ^^CO linewidths (Table 4) are typically 1-2 km s"^ and 
those of the C^^O lines, although narrower, range between 0.3 
and 1.6 km s"^ Several velocity components are present, either 
seen (in a few C^^O spectra) as two distinct peaks, or inferred 
from the non-Gaussian lineshapes, as shown by S6 spectra in 
Eig. 7. Assuming that the dust and gas in these cold dense clumps 
is in thermal balance, then the gas temperature is close to lOK 
and the sound velocity in CO is low, ~ 0.05 km s"^ . The molec- 
ular linewidths are therefore all suprathermal. It is interesting 
to compare the level of non-thermal support provided by these 
gas motions to the self-gravity. This can be done by compar- 
ing the gas linear mass density with its critical (or virial) value, 
mvir = 2cr^/G which is the largest mass per unit length of a self- 
gravitating cylinder, given the non-thermal support provided by 
the internal motions of dispersion cr = Av/2.35 (Eiege & Pudritz 
2000). 

The virial linear mass densities inferred from the C^^O 
linewidths of Table 4 range between lOMopc"^ for SIO, to 
26OM0pc"^ for SB. A crude comparison with the linear mass 
densities estimated at the resolution of the Planck data (Table 2) 
suggests that the ratios m/mvir are within a factor of a few of 
unity. Therefore, as a first approximation, these cold clumps are 
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Fig. 6. Spectral energy distribution of various components in 
the SI (upper panel) and S5 (lower panel) fields, including: the 
Planck SED according to the C3P0 catalogue (solid circles and 
thick solid line), aperture photometry with 27C and ?>1" or 58'' 
radius aperture sizes (squares and thin solid lines) and the back- 
ground SED (triangles and dashed line). The temperatures and 
P values derived from the modified blackbody fits are given in 
the figure (see also Table 3). The background SED is that of the 
median surface brightness in the reference annulus surrounding 
the larger aperture (see Fig. 5). 



in rough equilibrium between self-gravity and non-thermal sup- 
port. 

The column densities have been estimated from a single line- 
of- sight towards each object. When available (in the case of 
IRAM observations), we have combined the /=1^0 and 2^1 
transitions to derive the excitation temperature and the den- 
sity using a Large Velocity Gradient radiative transfer model 
(Goldreich & Kwan 1974). In the case of the Onsala and APEX 
observations, we only have data for one transition; in these cases, 
we assume an excitation temperature of 11 K and use the C^^O 
line intensity to estimate the column density of the cores. The 
molecular hydrogen column densities are calculated assuming a 
C^^O fractional abundance of [Ci^0]/[H2] = 10"^ . 

The results are presented in Table 4. The table includes esti- 
mates of the mean density, calculated assuming a core size cor- 
responding to a diameter of 74''. The column densities are close 
to the values derived earlier (Table 3) from the dust continuum 
data, using the smaller aperture sizes (i.e., 74" or 116"). The 
estimated average densities of the cores are only of the order 
of 10^ cm"^, which is lower than the densities found in nearby 
dense cores (e.g., Myers et al. 1983). The reason for this is that 
C^^O is not a good tracer of dense gas because of the significant 
degree of CO depletion onto dust grains. Within a cold core, 
the depletion could be almost complete (e.g., Pineda et al. 2010) 



properties of a sample of Galactic cold clumps 13 




I.. I p I ^ I N ^ 1 1 I I.I. I t ■ I I , ^ ■ I I .1 

-a -+ -2 -32 , -ED -Ifl 

LSR VWodtv Ckm T^) 

Fig. 7. CO spectral lines from the Planck clumps SI and S6 
(taken towards the central positions given in Table 1 , and with 
beam apertures in Table ??. The spectra are in the antenna tem- 
perature scale (T*). 

and so the investigation of this phenomenon towards the Planck 
cores is a natural topic for follow-up studies. 

5. Discussion 

5.1. The validity of the cold core detection method 

A perhaps unexpected result is the physical size of the Planck 
detections. They may have been anticipated to be compact point 
sources, the result of small cores being diluted within the large 
beam of Planck. However, they are instead found to be signifi- 
cantly extended and elongated, and embedded in filamentary (or 
cometary) larger-scale structures. As discussed in previous stud- 
ies, the boundary of a core/clump/cloud structure is always dif- 
ficult to assess. Curtis & Richer (2010) for instance have shown 
how much it depends on the method used to identify and extract 
the core parameters. 

In our case, the detection method used to extract sources 
from the Planck data is based on the colour signature of the 
objects, designed to enable us to detect the cold residuals af- 
ter the removal of the warmer background. This results in the 
discovery of a diff'erent, more extended cold component with 
a more complex morphology than sources found with methods 
that identify structures on the basis of surface brightness (e.g., 
the "clumpfind" algorithm of Williams et al. (1994), or a multi- 
scale wavelet analysis). The high sensitivity of Planck helps to 
better separate the warm and cold components, particularly in 
combination with data at somewhat shorter wavelengths. 

The faint, mostly filamentary, emission of cold dust that we 
detect here could perhaps be called a "cold matrix" linking the 
substructures to each other over a broad range of scales. By com- 
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Table 4. CO line parameters for selected lines-of- sight in seven 
cold core fields. 



-50 50 
R.A. offset (orcsec) 

Fig. 8. The S6 source spectral line intensity maps for C^^O 
/=1^0 (contour steps:0.5 Kkms"^ top image) and ^^CO 
/=1^0 (contour steps: 1 K km s"^ bottom image).The white 
disk indicates the 30-m HPBW size at the 2.7 mm wavelength, 
valid for both maps. 



bining Planck and Herschel it is possible to probe in detail the 
link between 10-pc scale filamentary structures and sub-parsec 
scale cold cores. Future studies with much larger samples should 
bring new insights to the origin of the cores and the mechanisms 
associated with their formation. 



5.2. The limitations of the SED fitting 

One has to be careful not to over-interpret SED fits based on a 
modified blackbody function, with a single colour temperature 
and optical depth. In reality the ISM has a broad distribution of 
dust temperatures and opacities (as well as column densities), 
within a Planck beam. Shetty et al. (2009) have discussed in de- 
tail the biases introduced by line-of- sight temperature variations 
and noise in the interpretation of such SED fits. Our study sug- 
gests that we are not badly contaminated by these eff'ects for sev- 
eral reasons. Firstly Planck's broad spectral coverage, combined 
with IRIS allows us to determining both the dust temperature 
and emissivity spectral index. And secondly, we stress that the 
cold core detection method identifies structures where the sub- 
millimetre SED is dominated by localised dust emission colder 
than its large scale background. 

An indirect indication that the parameters derived from the 
SEDs fits are not strongly biased is that the physical properties 
inferred cover the range of values known to be those of pre- 
stellar cores. Such cores have diff'erent properties whether they 
belong to regions of low- or high-mass star formation. Their col- 
umn densities range between a few times 10^^ or 10^^ cm"^ in 
nearby star forming regions (Motte et al. 1998; Kauff'mann et al. 
2008; Enoch et al. 2006, 2007; Hatchell et al. 2005; Curtis & 
Richer 2010; Andre et al. 2010) to ~ 10^^ cm'^ in the cores of 
infrared dark clouds (IRDCs, Simon et al. 2006; Rathborne et al. 



Target Telescope 


Line 


FWHM 








(kms-^ 


(cm ^) (cm ^) 



SI 



S3 
S6 

87^^ 

SB 
S9A 
S9B 
SIO 



IRAM 
IRAM 
IRAM 
IRAM 
IRAM 
IRAM 
Onsala 
Onsala 
Onsala 
IRAM 
IRAM 
IRAM 
IRAM 
IRAM 
IRAM 
IRAM 
IRAM 
Onsala 
Onsala 
Onsala 
Onsala 
Onsala 
Onsala 
APEX 
APEX 
APEX 



i^CO J=\- 
i^CO /=2- 
i^CO /=1- 
i^CO /-2- 
C^^O/^l- 
C^^O /=2- 
i^CO /=1- 
i^CO /=1- 
C^^O/^l- 
i^CO /=1- 
i^CO J=2- 
C^^O/^l- 
C^^O /=2- 
i^CO J=l- 
i^CO /-2- 

C^^O /=2- 
i^CO /=1- 
C^^O/^l- 
i^CO /=1- 
C^^O/^l- 
i^CO J=l- 
C^^O/^l- 
i^CO /=2- 
i^CO /=2- 
C^^O /=2- 



0.71 
0.78 
0.53 
0.50 
0.30 
0.32 
2.64 
1.71 
1.22 
1.44 
1.42 
0.82 
0.91 
1.60 
1.70 
1.40 
1.60 
2.06 
1.82 
2.54 
1.57 
1.95 
1.41 
1.39 
0.82 
0.35 



+ 0.04 
+ 0.07 
+ 0.01 
+ 0.02 
+ 0.02 
+ 0.01 
+ 0.03 
+ 0.04 
+ 0.08 
+ 0.02 
+ 0.02 
+ 0.01 
+ 0.01 
+ 0.10 
+ 0.10 
+ 0.20 
+ 0.20 
+ 0.07 
+ 1.07 
+ 0.04 
+ 0.12 
+ 0.05 
+ 0.11 
+ 0.02 
+ 0.03 
+ 0.20 



3000 3.0x 1021 

2700 4.0 X 10^1 

19000 6.7x 10^1 

5000 1.4x 10^2 

7300 5.6 X 10^1 

6400 1.3 X 10^2 

25000 1.8x 1022 

23000 1.6x 1022 

23000 1.8 x 1022 



^ Average density along the line of sight. 

^ Gas column density converted to H2 column density using the follow- 
ing conversion factors: A^H2 = IO^xNC^CO); and A^H2= lO^xA^(Ci^O). 

^ For this source, the lines are an assembly of three components, with 
strongly asymmetric shapes, which are difficult to analyse separately. 



2010; Peretto et al. 2010). These IRDCs, discovered by means of 
mid-infrared absorption towards the bright background emission 
of the Galactic Plane (MSX and ISOGAL surveys, Egan et al. 
1998; Perault et al. 1996), are thought to be the sites of formation 
of massive stars and star clusters (Rathborne et al. 2006). The 
inferred average densities depend on their size, which ranges be- 
tween < 0.1 pc for low-mass dense cores of ~ 1 M©, to 0.5 pc 
or more for -10^ M© IRDC cores. Our Planck cold clumps sub- 
sample is not located within the MSX IRDCs spatial distribution, 
so we cannot compare them directly with any IRDC association. 
However, the properties we derive fall well within this range. It is 
noteworthy that the Planck cold clumps are elongated structures 
which tend to exist within extended filaments (up to 30 pc in 
the present subsample), with thickness up to ~ 3 pc. The IRDCs 
identified as extinction peaks in mid-infrared maps have similar 
lengthscales and column densities, but they are usually thinner 
and therefore have higher densities. It is also interesting that the 
linear mass densities of the filaments in which the Planck cold 
clumps are embedded cover the same range as the warm filamen- 
tary structures associated with active star-forming regions, i.e., 
extending up to several 100 Mopc"^ 

The temperatures derived in the sample studied here cover 
a range (10-14.5 K) similar to previous estimates on a num- 
ber of cold condensations detected using multi-wavelength 
submillimetre observations with the balloon-borne experi- 
ments: PRONAOS, r ~ 12 K in star-forming and cirrus regions 
(Stepnik et al. 2003; Dupac et al. 2003; Bernard et al. 1999); 
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ARCHEOPS, 7-18 K (Desert et al. 2008); and BLAST, 9-14 K 
(Netterfield et al. 2009). These temperatures are all lower, how- 
ever, than those measured in the so-called "quiescent" IRDC 
cores (Rathborne et al. 2010) that span from 17 to 30 K for a 
similar range in column density (0.3-3x10^^ cm"^) and mass 
(10-10^ Mq). The Planck cold clump population may therefore 
be representative of a still earlier stage of evolution of cold dense 
cores. 

6. Summary and perspectives 

We have presented a preliminary analysis of a sample of 10 
sources from the C3P0 catalogue in order to illustrate and bet- 
ter probe the nature and properties of the cold objects detected 
with Planck. The sources have been chosen to span a broad 
range in temperature, density, mass and morphology (inluding 
filaments and isolated/clustered structures) in a variety of envi- 
ronments, from star-forming regions (both remote and nearby), 
to high Galactic latitude cirrus clouds. The main findings are as 
follows: 

1. The sources are significantly larger than the Planck beam, 
with elongated shapes, and appear to belong to filamentary 
structures, with lengths up to 20 pc; 

2. The physical parameters of the sources have been de- 
rived from SEDs by combining Planck (HFI bands at 
857, 545, and 353 GHz) and IRIS data, finding T-IO- 
15 K (with a mean value of 12.4 K), 1.8-2.5 (mean 
2.2) and A^H2 ~ 0.8-16 x 10^^ cm"^, from which we infer 
linear mass densities in the range m = 1 5-400 M0pc"^ 
masses M~ 3.5-1 800 M©, bolometric luminosities L~l- 
300 Lo, and L/M ~ 0.1-0.9 Lq/Mq; 

3. Except for the faintest source (which lies at high Galactic 
latitude), a clear signature of cold dust emission is visible di- 
rectly in the 1° X 1° maps of dust temperature, spectral index 
and column density, with colour temperatures of typically 
~ 14 K, surrounded by a warmer extended (often elongated) 
emission at around 16-18 K; 

4. Herschel and AKARI observations at higher angular resolu- 
tion have revealed a rich and complex substructure within the 
Planck clumps, in most cases the substructures being colder 
(down to 7 K) than the Planck-dctcctcd clumps, although in 
two cases, the substructures are warmer because they har- 
bour compact infrared objects, likely protostellar sources at 
an early stage; 

5. Molecular line observations of 7 of the sources show that 
all of them are clearly detected in ^^CO and C^^O, ^^CO 
linewidths typically l-2kms"^ and C^^O lines always nar- 
rower (down to ~0.3kms"^) but still clearly suprathermal, 
given the anticipated low temperature of the gas, suggesting 
that the support of Planck cold clumps against self-gravity is 
dominated by non-thermal motions. 

Although we have focussed here on a very small sample of 
10 clumps, the results are indicative of what might be expected 
from the more ambitious studies which will follow. We have 
already shown that the C3P0 list of Planck cold sources con- 
tains objects with a wide variety of physical properties. These 
are probably associated with diff'erent evolutionary stages of the 
star formation process, from quiescent, cold and starless clumps, 
through prestellar stages to very young protostellar objects still 
embedded in their cold surrounding cloud. Forthcoming pa- 
pers will present a more detailed analysis of the sources using 
Herschel data combined with Planck for larger and more stat- 
tistically robust samples. More detailed physical modelling will 



allow for characterisation of clustering and fragmentation within 
the Planck clumps, addressing the question of their evolutionary 
stage, along with the study of stability of the starless substruc- 
tures. 

The first unbiased all- sky catalogue of cold objects provided 
by Planck off'ers the opportunity to investigate the properties of 
the population of Galactic cold objects over the entire sky. The 
full sample will include objects at the very early stages of evolu- 
tion, in a variety of large-scale environments, and in particular, 
outside the well-known molecular complexes, at high latitude 
or at large distances within the Galactic Plane. For this purpose, 
dedicated follow-up observations are needed in both higher reso- 
lution continuum mode and spectroscopy, which is the objective 
of the "Galactic Cold Cores" key programme, planning a follow- 
up with PACS and SPIRE of about 150 Planck cold clumps. In 
parallel, similar and/or complementary follow-up studies will be 
possible on the basis of the Early Cold Core catalogue which 
has been delivered to the astronomical community, providing a 
robust sub-sample of the C3P0 catalogue, with more than 900 
Planck cold clumps distributed over the whole sky. 
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Appendix A: Selected sources, illustration of 
detection method, source elongation and 
observational parameters 

We list below the main features of the selected fields as known 
based on previous studies. 

51 lies in a tenuous high-latitude cloud, MCLD 126.5+24.5, 
located at the border of the Polaris Flare, a large molecular cir- 
rus cloud in the direction of the north celestial pole (Heithausen 
& Thaddeus 1990), at an estimated distance of 150 pc (Bensch 
et al. 2003). The cloud was first noted from the POSS prints 
by Lynds (1965) as the small reflection nebula LBN 628. Its 
cometary globule shape appears similar to what is usually found 
in active star formation regions, although this nebula is far from 
any such region. Boden & Heithausen (1993) have studied the 
gas properties of the densest part of the cloud, observing transi- 
tion lines of the dense tracors ^^CO, ^^CO, H2CO, and NH3. The 
ammonia core (~0.14pc x 0.07 pc) corresponds to 0.2Mo and 
a density of 4000cm"^. Although this density is high compared 
to the cirrus-like environment, Boden & Heithausen (1993) find 
that both the cloud and the core are probably gravitationally un- 
bound. The line ratios found are consistent with shocked gas be- 
ing compressed by the North Celestial Pole Hi loop (Kramer 
et al. 1998). 

52 is embedded in the tail (~ 30' south) of the well- studied 
cometary globule CGI 2. This isolated globule is located at high 
Galactic latitude (/ = 316.5°, = 21.2°), with a 1° long nebu- 
lar tail nearly perpendicular to the Galactic Plane. Its distance 
of 550 pc has been determined using a photometry-extinction 
method by Maheswar et al. (2004). Molecular studies using 
i^CO (van Tin et al. 1975; White 1993; Yonekura et al. 1999), 
C^^O (Haikala et al. 2006), H2CO (Goss et al. 1980) and 
NH3(Bourke et al. 1995) reveal star formation activity, with the 
presence of a highly collimated bipolar outflow in the head of the 
cloud. CGI 2 seems to be a rare example of triggered star forma- 
tion at relatively large Galactic height (Maheswar et al. 2004). 

53 is located in the constellation of Cassiopeia and has no 
counterpart in the SIMBAD database. On the sky, the closest 
Lynds catalogue sources are LDN 1358 and LDN 1355, which 
are, respectively, at distances of 114 and 118'. These Lynds 
sources are associated with the Cepheus flare at a distance of 
200+50 pc (Obayashi et al. 1998; Kaufl'mann et al. 2008). The 
Planck detection is associated with an AKARI FIS bright source 
catalogue source that has a cirrus type spectrum (Yamamura 
2008). 

54 is in the Orion region, within a dark cloud that was 
mapped by Dobashi et al. (2005) in extinction, derived from 
Digitized Sky Survey images. The source is located six degrees 
east of the M42 nebula, in the Dobashi et al. cloud 1490, close 
to its clump P43. It has not been the subject of any dedicated 
studies, so far. In higher resolution extinction images derived 
from the colour excess of 2MASS stars, cloud 9711 is seen to 
reside within a narrow filament. In the Planck data, the filament 
is clearly visible, but remains unresolved. The source itself is a 
compact, slightly elongated clump that stands out just as well 
in the individual Planck frequency maps as in the cold residual 
map. The Nanten CO data show a peak at the same position, with 
a line intensity close to 20 K km s"^ 
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55 is a filamentary structure located near the Vela Molecular 
Ridge (hereafter VMR), a giant molecular cloud complex in 
the outer Galaxy. Murphy & May (1991) mapped the VMR 
in CO, but the mapping did not cover S5. However, Otrupcek 
et al. (2000) carried out CO observations of the dark cloud DCld 
276.9+01.7, with angular extent of 16' by 3', and this is inter- 
posed on the filament. The detected CO line is rather strong 
(r^ = 6.1 K), with a velocity (relative to the local standard 
of rest) of ~2kms"^. We derive a kinematic distance for the 
cloud of ~ 2 kpc. This is consistent with the distance estimate 
of the VMR: Liseau et al. (1992) give a photometric distance of 
(0.7 ± 0.2) kpc for parts A, C and D of the VMR and 2 kpc for 
part B, using the nomenclature established by Murphy & May 
(1991) and agreeing with their results. Thus we consider it is 
likely that Planck-4110 is part of the VMR. Several studies have 
been published on the VMR young star-forming region, for ex- 
ample, Olmi et al. (2009) found 141 BLAST cores (starless and 
proto- stellar) in the Vela-D region. 

56 has been detected in the anti-centre direction (/, 
Z7=176.18, b = -2.11), at a distance estimated to 2 kpc using 
the extinction signature method by Marshall et al. (2006). The 
filament- shaped core seen with Planck does not coincide with 
any known object, but is likely associated with the Perseus arm. 

57 has been observed by Ungerechts & Thaddeus (1987) as 
part of their large CO survey of Perseus, Taurus and Auriga. 
It is included in their 12th area (see their Table 1), covering 
41.8 deg^, with a distance of 350 pc. For this large region, they 
estimate a virial mass that is four times higher than the CO 
mass, showing that this region is out of equilibrium. The Planck 
source is not associated with any known objects in the SIMBAD 
database.. 

S8, S9 and SIO correspond with Herschel SDP sources 
PCC288, PCC550 and PCC249, respectively. These fields were 
selected in September 2009 from the cold core detections in 
the First Light Survey of Planck, in order to perform the first 
follow-up observations of our Herschel key-program during the 
Herschel SDP. The fields were mapped with both PACS and 
SPIRE (with map sizes from 18' to 50'), and the observations 
have been described in detail in dedicated papers (Juvela et al. 
2010, 2011). 

Field S8 is located in Cepheus, at the interface between the 
Cepheus F molecular cloud (Sargent 1977) and the young stel- 
lar group Cep 0B3b (Jordi et al. 1996), which is one of the 
youngest nearby stellar groups. It has been suggestsed that this 
region could show direct evidence for star formation triggered 
by the OB association. S9 is an active star formation region 
close the S140 Hii region. Planck has provided two detections, 
PCC288-P1 and PCC288-P2 (see Juvela et al. 2010), with the 
southern source PCC288-P2 being the colder one. The two cores 
have also been identified in CS /=1^0 mapping of the area by 
Tafalla et al. (1993). The field SIO is part of the Musca cloud fil- 
ament, where Planck shows two secure detections of cold cores 
(r ~ 1 1 K). The cores have been studied previously in ^^CO and 
C^^O by Vilas-Boas et al. (1994). 
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Fig.A.5. Comparison of the source elongation with the beam 
shape determined with the "FEBeCoP" tool (Mitra et al. 2010). 
The ellipses correspond to the estimated size of the Planck cold 
clumps, and the dashed circles trace the local beam shape. 
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Fig. A.l. Multi-band emission maps of the sources SI, S2, S3, and S4, 
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Fig. A.2. Multi-band emission maps of the sources S5, S6, S7, and S8. 
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Fig. A.4. Illustration of the steps used in the detection and core extraction methods, as applied to source SI. The first column shows 
the initial surface brightness emission maps in the IRIS and HFI bands. The rows show the process applied to the 100 yum, 857, 545, 
and 353 GHz data, respectively. As described in Paper I, the method follows the following main steps: a "Warm Ref" is built from 
the IRIS lOOyum map and extrapolated to the Planck bands using the local colour estimated around each pixel of the Planck map 
(these colour maps are shown in the second column); the 'Warm Bkg' map obtained at a given frequency is then removed from the 
Planck map, revealing the "Cold Residual"; and finally an elliptical Gaussian fit is performed on the colour maps to derive core the 
size parameters, while aperture photometry is performed on the HFI bands, integrating the signal within the elliptical Gaussian. 



